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We generalize our previous work on the phase stability and hydrodynamic of polar liquid crystals 
possessing local uniaxial Coo^-symmetry to biaxial systems exhibiting local C2„-symmetry. Our 
work is motivated by the recently discovered examples of thermotropic biaxial nematic liquid crystals 
comprising bent-core mesogens, whose molecular structure is characterized by a non-polar body axis 
(n) as well as a polar axis (p) along the bisector of the bent mesogenic core which is coincident with a 
large, transverse dipole moment. The free energy for this system differs from that of biaxial nematic 
liquid crystals in that it contains terms violating the p — » p symmetry. We show that, in spite 
of a general splay instability associated with these parity-odd terms, a uniform polarized biaxial 
state can be stable in a range of parameters. We then derive the hydrodynamic equations of the 
system, via the Poisson-bracket formalism, in the polarized state and comment on the structure of 
the corresponding linear hydrodynamic modes. In our Poisson-bracket derivation, we also compute 
the flow-alignment parameters along the three symmetry axes in terms of microscopic parameters 
associated with the molecular geometry of the constituent biaxial mesogens. 

PACS numbers: 61.30.-v, 47.20.-k, 03.50.-z 



I. INTRODUCTION 



In the past few decades, the field of liquid crystallog- 
raphy has benefited from a tremendous progress in the 
increasing sophistication of technology and experimental 
designs. The ability to manipulate fine microscopic de- 
tails of constituent molecules, such as geometrical shape, 
electric and magnetic moments, as well as chemical affin- 
ity, has enabled the creation of many interesting liquid- 
crystalline phases of ever increasing complexity, both in 
terms of their characterization and properties. Since 
Freiser's seminal theoretical work [1\ predicting that com- 
pounds composed of biaxial molecules should in principle 
engender a biaxial nematic phase, and the experimental 
discovery of the lyotropic biaxial nematic phase a decade 
later by Yu and Saupe , the experimental realization of 
a thermotropic biaxial nematic phase has long been the 
holy grail of liquid crystallography. Aside from the poten- 
tial technological applications, biaxial nematic systems 
also hold theoretical interest in that they provide the sim- 
plest example of an ordered medium with a non-Abelian 
fundamental homotopy group [31. In the ensuing decades, 
there were several claims [j, la, @ to the discovery of 
systems displaying biaxial nematic phase as confirmed 
by optical techniques. However, under the scrutiny of 
the more reliable method of NMR spectroscopy Q , these 
prior claims were proven to be false hopes. In the mean 
time, novel engineering attempts were made involving 
more complex systems, such as binary mixtures of hard 
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rods and plates, that exhibit a macroscopic phase with 
two symmetry axes [8j, where each species in these bi- 
nary mixtures orientate in different directions that are 
orthogonal to each other. Finally, in 2004 unequivocal 
evidence materialized demonstrating the existence of the 
thermotropic biaxial nematic phase based on the tech- 
niques of NMR spectroscopy [3, [l(| and X-ray scatter- 
ing In these experiments, the elusive biaxial nematic 
phase was realized by systems composed of bent-core 
mesogenic molecules, which had first garnered attention 
for the discovery by Niori et al. [12] because they exhib- 
ited unusual smectic phases with ferroelectric and antifer- 
roelectric properties. It has been postulated [1, [l(| that 
strong intermolecular associations that originate from the 
large electric dipole moment pointing towards the center 
of these "V" -shaped mesogens, rather than their pure ge- 
ometry, should play a crucial role in the stability of the bi- 
axial nematic phase. One strong evidence supporting this 
claim comes from the experimentally observed apex an- 
gle, formed by the two arms of each mesogenic molecule, 
of 140 degrees which is far from t he p redicted theoretical 
optimum value of 109.4 degrees Thus, it is an 

important consideration that one of the two molecular 
symmetry axes of these bent-core mesogens, along which 
the direction of the dipole moment lies, is polar. 

Such general issues on the geometry of mesogenic 
molecules and the stability of the resulting macroscopic 
phases have been investigated by a few numerical stud- 
ies [l3|, [3 [34], [62j • Specifically, it has been found that 
when the angle between the arms of these bent-core 
molecules becomes the tetrahedral (which has generally 
not been observed in mesogenic systems exhibiting the 
biaxial nematic phase), the isotropic phase undergoes 
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FIG. 1: A schematic representation of the bent-core meso- 
ecn with orthonormal molecular axes {£>i, £2, £3}. The 
molecular dipole points along Ox, while the molecular body 
axis lies along O3; O2 points into the page. Macroscopically, 
we define the orthonormal triad {n, p, m} with the polar and 
nematic director fields along p and n, respectively. 



a continuous phase transition into the biaxial nematic 
phase [HG3. The effect of the anisotropy in the length 
of the two arms as well as the angle between them has 
been investigated in a more recent Monte Carlo simu- 
lation [11}. Analytically, a few other theoretical issues 
have also been looked at, including static properties such 
as stability of phases fl6L [Til WsL fl9l [20 . l2lj ] , elastic- 
ity and viscosity [H, l23l . \24 . 12511 , and critical behav- 
ior [26|. The dynamical properties [27|, [H, [U, M, SH [H 
of biaxial nematic systems have also been considered. 
For the most part, many of the aforementioned the- 
oretical works have been mainly focused on biaxial 
systems exhibitingthe macroscopically non pol ar D^hr 
symmetry Q, [H, M, M M M M M M M HI , in 
which both of the symmetry axes are non-polar. 

In this paper, we focus on the stability and dynam- 
ics of the polar biaxial phase (C*j„) that preserves the 
symmetry of the molecular polar structure of these bent- 
core mesogens. The alternative scenario of macroscopic 
polarity in a biaxial nematic phase, wherein two symme- 
try axes are polar, has been considered in plj . To fully 
characterize this bipolar biaxial phase, the spatial orien- 
tation of the molecules must be specified by the full set 
of Euler angles, thereby generalizing the corresponding 
distribution function from the spherical harmonics {l 7 ^} 
to the Wigner D-functions { D^ n m , } [35] , which form a 
complete set of orthogonal functions with respect to in- 
tegration over SO(3). The necessary order parameters, 
of their bipolar biaxial nematic phase, that exhibits local 
C s -symmetry, then reduce naturally to the components 
of two vectors. In more general instances, the most di- 



rect consequence of biaxiality is the generalization of the 
alignment tensor to contain the usual Maier-Saupe order- 
parameter as well as a new biaxial order-parameter. For 
concreteness, we define our coordinate systems as shown 
in Fig. Q] We assume that the dipole moment and the 
body axis of a bent-core mesogen point along its molec- 
ular axes i>iand £3, respectively. The third unit vector 
identifying the molecular coordinate system is £2= 

I>1, SUCh that J^a "a^aj = 

For general systems exhibiting biaxial symmetry, achi- 
ral ordered phases are fully characterized by both the 
second-rank biaxial alignment tensor and the align- 
ment vector P. These in turn can be written as 



,(r) = S(r) 



nj(r)n,-(r) 
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+T(r) \pi(r)pj(r) - m^jm^r)] , (1) 

P(r) = S p (r)p(r)+S„(r)n(r)+S m (r)m(r) , (2) 

where S = (P 2 (cos6»)} and T = (sin 2 6> cos 20) , with 
and <f> denoting the polar and azimuthal angles, respec- 
tively. The unit vectors {p, m, n} are an orthonormal 
triad. In [f5j], Lubensky and Radzihovsky introduced a 
third-order tensor parameter for the bent-core molecules 
in order to fully characterize the chiral phases also ob- 
served in these systems. Since we do not consider such 
chiral phases in our work, it suffices to consider only the 
alignment tensor and vector as defined in Eqs. ([I]) and 
©■ 

In the polar biaxial phase where the symmetry of the 
macroscopic phase coincides with that of the constituent 
molecules (C^u), we assume that, on average, the non- 
polar body axis of the mesogens points along n(r), while 
the polar axis induced by the microscopic dipole moments 
lies along the unit vector p(r). Thus, S n = S m = 0. Fur- 
thermore, n(r) and p(r) are always orthogonal by con- 
struction and m = n x p. The vector m is therefore not 
an independent quantity. The biaxial alignment tensor 
can be rewritten solely in terms of p and n as 



= (S + T) ( TliTlj - ^Sij 



2T p iPj - -6ij 



(3) 



The scalar quantity S is the usual Maier-Saupe order 
parameter that describes alignment along n(r), the di- 
rection corresponding to the maximum eigenvalue of 
UiQijUj = ^S. The biaxial order parameter T describes 
ordering of the mesogens in the m — p plane. The di- 
rection along p corresponds to the second-largest eigen- 
value of PiQijPj — —^S + T. The scalar field S p en- 
codes information on the degree of local alignment of 
the dipoles [361 ] . The characterization of biaxal nematic 
phases with different local symmetries would generally 
require different sets of order parameters. 

A comprehensive classification of possible phases 
for bent-core molecules has been given by Luben- 
sky and Radzihovsky [lj| and also by Mettout [2fj| . 
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With the identification of the proper order param- 
eters, group-theoretic arguments and Landau mean- 
field theory can be used to catalogue the universal- 
ity classes of these phase transitions. It has been 
found in [151 ] that transitions to the polar biaxial phase 
with C2i,-symmetry {S p ^ 0, S ^ 0, T ^ 0} are allowed 
from various higher-symmetry phases, including the po- 
lar unaxial (Coo„) phase {S p ^ 0, S ^ 0, T = 0} and the 
uniaxial nematic (-Doo/i) phase {S p = 0, S ^ 0, T = 0}, 
as well as the nonpolar biaxial nematic (I?2h) phase 
{S p = 0, S ^ 0, T 0}. Assuming that the nematic di- 
rector field aligns along n and the polar director field 
along p, one can envision the transition to the polar bi- 
axial (C l 2 t ,)-phase from the polar uniaxial phase by the 
development of biaxial order in the plane perpendicu- 
lar to p, or similarly, from the uniaxial nematic phase 
via the establishment of a vector order along p. Rota- 
tional invariance in the plane normal to the director n 
implies that both of these transitions will be in the XY- 
universality class [H, H3- O n the other hand, the tran- 
sition to the polar biaxial (C^u) phase from the nonpolar 
biaxial nematic (£*2/i) phase is characterized by the de- 
velopment of vector order along one of the three twofold 
symmetries already specified in the biaxial nematic phase 
and thus renders the transition a member of the Ising uni- 
versality class with well-known critical properties 15 , [37| • 

In what follows, we will investigate the stability and 
the dynamics of the polar biaxial phase (C% v ), in the 
linear hydrodynamics regime, thereby generalizing the 
results of our prior work [36| on the uniaxial case. We 
derive the hydrodynamic equations from a microscopic 
model using the Poisson-bracket formalism and obtain 
expression for the flow- alignment parameters. The struc- 
ture of the hydrodynamic equations in the biaxial case 
differs qualitatively from the unaxial case, since biaxial- 
ity breaks completely the rotational group SO (3). This 
symmetry breaking generates a total of three Goldstone 
modes, in contrast to only two found in the uniaxial case 
due to the residual SO(2) group. The structure of the 
corresponding hydrodynamic modes is also more com- 
plicated due to the discrete orthorhombic symmetry in 
the biaxial mesogens, which is a lower symmetry group 
than the cylindrical symmetry found in uniaxial meso- 
gens. The existence of polarity again implies additional 
parity-odd terms in the free energy that would otherwise 
be forbidden in the nonpolar case. The most important 
of such terms couple density fluctuations and polariza- 
tion splay and render the uniform polarized biaxial phase 
unstable in a range of parameters. 

Our derivation starts in section II where we introduce 
a microscopic definition of the various order parameters. 
In section III, we consider the stability of the polarized 
phase from the viewpoint of the free energy. In par- 
ticular, we show that polarity can be compatible with 
orientational order in liquid phases. We then obtain the 
hydrodynamic equations for the biaxial polarized phase 
in section IV. We consider these equations in the linear 
hydrodynamic regime and comment on the structure of 



the resulting hydrodynamic modes in section V. Section 
VI concludes our paper. 



II. DEFINITIONS OF MICROSCOPIC FIELD 
VARIABLES 

In our previous work [36|, we studied the hydrody- 
namics of polar liquid crystals with local Coot, -symmetry. 
We applied the Poisson-bracket formalism [H, [39[ to 
derive the hydrodynamic equations for the conserved 
variables — mass and momentum densities — as well as for 
the polar director field. We obtained a microscopic ex- 
pression for the flow alignment parameter A [13, S3, EH, 
|42| and showed that the hydrodynamic equation for the 
polar director field have the same form as its nematic 
counterpart, albeit with a different A [H, 0, [H, HU . As 
mentioned in the Introduction, we have introduced the 
polarization P as the order parameter characterizing the 
polar liquid phase in uniaxial systems. In the case of a 
biaxial phase in which one of the two director fields is 
polar, we will need to employ the alignment tensor Qij 
in conjunction to the alignment or polarization vector P 
in order to fully characterize the polar phase. 

We consider a system composed of N identical biaxial 
mesogens, indexed by a. Based on the three-atom model 
in |l5( | . we assume that each molecule consists of three 
atoms, treated as point particles, along the two arms of 
the boomerang, as shown in Fig. (UJ). Each molecule is 
described by a molecular coordinate system {£>i, j> 2? £3} 
as defined in Fig. (fT]). The position of the center of mass 
of the a-th molecule is given by 



M 



2m + M 



(4) 



The atoms of mass m have charges q, and the central 
atom of mass M possesses charge Q. We assume that 
molecules are neutral, 



2q + Q = Q 



(5) 



but possess a no n- vanishing dipole moment d a , given by 



d a = QL a cos^ Q 0" 



(6) 



where is the unit vector along the polar molecular 
axis (Fig. [JJ . We note that while the molecular dipolar 
interactions are physically important in stabilizing the 
observed biaxial liquid-crystalline phase, the asymmetry 
in the molecular shapes alone also gives rise to polarity in 
the system. The alignment vector P can then be defined 
as the coarse-grained sum of the microscopic polar axes, 



P r = rr 



/>?<5(r-r") 



(7) 



where mo = M + 2m, p(r) is the local mass density, and 
can be expressed in terms of a unit vector p, as given in 
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Eq. {2J with, S n — S m — 0. Analogously, the alignment 
tensor is defined in terms of a coarse-grained sum of the 
traceless, second molecular mass-moment tensor Rf.\ 



R(r) 



££g<5(r-r<*) 



where 



TDOi 



E 1 



ArfAr?" - (Ar^) 



(8) 



(9) 



Here fi labels the atoms within each molecule and 
Ar"' 1 = r aAI — Ar" is the position of each atom relative 
to the center of mass. The molecular mass distribution 
tensor R\j can also be written with respect to the basis 
of the orthonormal molecular coordinate system as 



£>a \ nana 



(10) 



as well as the anisotropy AI a 



AI a = 1° , 



2m (L a f sin 2 ^ a 



mM 
2m + M 



(L a ) cos 2 V c 



(20) 



Subsequently, we neglect fluctuations in the molecular 
shape and thus assume that R" is thus independent of 
a. We can then write 



R(r) 



AI 

m ' 



p(r) 



(21) 



In the limit of if) — > 7r/2, its molecules become needle-like, 
and we recover the expected ratio of I j AI = 2 [43l l6l|. 

The biaxial alignment tensor in Eq. ([5]) can thus be 
identified with the macroscopic expression in Eq. ([1]) 
expressed in terms of orthogonal unit vectors m, n, and 



with 



= 1,2,3, (11) 



and 



R% 



-R% 



mM .(L«) 2 cos 2 ^, 



2m + M 



(12) 



2m {L a ) 2 sin 2 ^ - 



mM (L°) 2 cob 2 ^. 
2m + M ' 



Finally, the scalar field R(r) is given by 
R(r) = J2KS(r-r a ) 

a, a 



(13) 

(14) 
(15) 



Equivalently, we can express the coefficients {R"} m 
terms of the molecular moment-of-inertia tensor 

Zg = m ^ [(Ar tt ^) 2 5 i:j - Ar^ArH , (16) 



whose components along the principal axes are related 
to the second molecular mass-moment coefficients as fol- 
lows: 



ret 

na i 33 



-BS 



T a _ J 33 _ ja _ Ji 33 
hi 2 ~h2 2 



(17) 
(18) 



We can then write down, in a straightforward manner, 
the trace I a 

I a ee Tr/g = Jf 1 +J 2 Q 2 +J 3 Q 3 

= 4m (L a f sin 2 ^ + „ 4mA * (L a f cos 2 ^ , 



2m + M 



(19) 



III. FREE ENERGY AND STABILITY OF THE 
POLARIZED PHASE 

Before proceeding with the derivation of the hydro- 
dynamic equations for the polarized state, it is instruc- 
tive to discuss the coarse-grained free energy of a polar, 
biaxial fluid. Compared with the polar unaxial phase, 
the free energy becomes much more complex due to the 
overall lower symmetry of the biaxial phase. We will 
again divide the free energy into a parity-even contri- 
bution that respects the symmetries of n — > — n and 
p — > — p and a parity-odd contribution that arises from 
the polar nature of p, F = F cvcn + F odd . The num- 
ber of distinct elastic constants is dictated by the dis- 
crete symmetries of the system. Starting with a general 
elastic free energy with terms quadratic in the gradients, 
i.e. terms of the forms (ViUj) (Vin m ), (ViPj) (V;p m ), 
and (Vjfij) (V;p m ), and minding symmetry considera- 
tions, we find that there are in total fifteen indepen- 
dent terms in the elastic free energy of biaxial nematics. 
Using the anholonomity (Mermin-Ho) relations for the 
three broken rotational symmetries [47|, |48| , which first 
appeared in the context of superfluid phases of helium- 
3 [4{J [5(| [U, [H, HH , three of the fifteen terms can be 
shown to be surface contributions [27| ■ The Mermin- 
Ho relations, which reflect the consequence of the non- 
commutative nature of finite rotations, can be expressed 
in terms of their locally-defined infinitesimal generators 
5® = (SQ 1 ,Se 2 , SQ 3 ) = (m • Sn, - (p • Sn) ,p ■ 5m): 



[8 1 S a -S a Si)Q = (<5i0) x (6 2 &) 



(22) 



where Si and 82 represent any first-order differential op- 
erators, such as dt or V. For example, to the first-order 
approximation of an infinitesimal variation of the direc- 
tor n, defined with respect to the free-energy functional, 
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The remaining twelve terms in the elastic free energy mass terms that are invariant under parity are given by 
are bulk contributions. The bulk elastic terms and the 



^° von = \ J d 3 r jd (^j + C 2 (Vp) 2 + K nl (V • n) 2 + K n2 (n • V x n) 2 + K n3 (n x V x n) 2 + K pl (V • p) 2 

+K p2 (p • V x p) 2 + +K p3 (p x V x p) 2 + K ml [V • (p x n)] 2 + K m2 [(p x n) • V x (p x n)] 2 
+K m3 [(p x n) x V x (p x n)] 2 + K mp [(p x n) • (V x p)] 2 + K pn [p • (V x n)] 2 

+K nm [n-(V x (pxn))] 2 } , (23) 

I 



where dp = p — po is the local fluctuation in density. 
The elastic free energy can be written down in many 
equivalent forms. The first nine elastic constants, from 
K n \ up to K m3 , have the simple interpretation of de- 
scribing the bulk deformations (splay, twist, and bend) 
along the three symmetry axes, respectively. The last 
three terms have no simple analogues in uniaxial elastic- 
ity. Theoretical calculations based on the microscopic ap- 



proach [22|, HH HH have been done to determine the mag- 
nitudes of these elastic constants from first principles, as 
the results strongly suggested a match to what one would 
intuitively expect that these "cross" -terms should take on 
smaller magnitudes than the other nine elastic constants. 

The existence of polarity allows for terms that break 
the p — > — p in the elastic free energy. The lowest-order 
contribution to the parity-odd free energy has the form: 



F° 



(id 



r/V<!D^(V-p) + D p (V-p) 

Po 



(V-n) 2 + 0(V 4 ) 



(24) 



By construction, Eq. (|24p involves only odd orders of the 
polar director p. As in the unaxial case, terms ~ V-p(r) 
give rise to surface terms that would favor a splay distor- 
tion of the mesogens, though the instability can be sup- 
pressed assuming appropriate surface stabilization. In 
addition to possible surface instability, it is apparent from 
the form of the new parity-odd i?-term that a bulk insta- 



bility can also develop. This intrinsic bulk instability to 
splay deformations has been well known in the uniaxial 
case [361 ] . We can understand qualitatively the existence 
of this bulk instability by rewriting the terms in the free 
energy involving density fluctuations and splay deforma- 
tions of the director field: 



(25) 



where <5p = p po is the fluctuation of the polar di- 
rector field about its equilibrium orientation. Thus, the 
existence of polarity along a symmetry axis in the biaxial 
phase, which allows for new parity-odd terms in the free 
energy, yields an overall downward renormalization of the 
splay elastic constant K p i . We emphasize again that the 
molecular asymmetry of the bent-core molecules is suf- 



ficient to induce the aforementioned polarity. Though 
in actual experimental realization of the biaxial phase 
(S p = 0,S'^0,T^0), the existence of a strong molecu- 
lar dipole moment and the resulting dipolar interactions 
are conjectured to be the essential contributing factor to 
the stability of this phase. Our discussion remains valid 
as long as polarity exists in the system, regardless of the 
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nature of its origin. 

While completing the square in the free energy expres- 
sion of Eq. (f25f demonstrates that the splay instability 
is a general bulk property of the polarized state, this sim- 
ple argument does not capture the anisotropic nature of 
the instability that will become apparent below from the 
analysis of the hydrodynamic modes of a polar biaxial 
system to be carried out in section V. The anisotropy of 
the splay instability is also fully apparent form a purely 
static calculation of the density-density correlation func- 
tion of the system. Such a calculation must be carried out 
enforcing the orthonormality constraint that n • p = 
which is not incorporated in the free energy. 



The free-energy as written in Eq. (|25|) indicates that 
splay fluctuations become energetically favorable and 
destabilize the ordered state when \B\ > K V \C\. The 
instability is suppressed if the fluid is incompressible 
since the renormalization of the Frank constant K p \ 
vanishes when C\ — > oo. The stable ground state for 
1-5 1 > \]K p \C\ is expected to be characterized by direc- 
tor configurations that are spatially inhomogenous in a 
direction that is perpendicular to the polar, broken sym- 
metry director p(r), with associated spatial structures in 
the density 

We can examine the possible relevance of the splay 
instability to experimental systems by performing nu- 
merical estimates of the various elastic constants. We 
consider a fluid of mesogens of linear size / interact via 
a short-range repulsive interaction of strength U. The 
mesogens forming the biaxial phase possess a permanent 
dipole moment of magnitude d that yields a dipolar cou- 
pling of strength Ud ~ k c d 2 /r 3 between two mesogenic 
molecules at distance r, where k c is the Coulomb's con- 
stant. The splay Frank coefficient has dimensions of en- 
ergy/length and can be estimated as K pl ~ U/l, and 
the bulk compressional modulus C\ has dimensions of 
energy density. Its magnitude is controlled by the repul- 
sive part of the interaction and thus can be estimated as 
C\ ~ U/l 3 . The elastic constant B is controlled by inter- 
actions that favor alignment of the mesogenic molecules 
and would otherwise vanish if the molecules are nonpolar. 
Thus, we estimate B ~ Ud/1 2 , where Ud is the strength 
of the dipolar interaction. In our estimate, splay insta- 

— 1/3 

bility will occur when Ud ~ U . Since r ~ p , we 
take po ~ 10~ 27 m and d ~ 4 Debye 0, [ll[ and estimate 
Ud ~ 0.01 eV. While there are no reported measure- 
ment of the elastic constants for thermotropic biaxial liq- 
uid crystals, theoretical analysis in [55| has established a 
mapping between the elastic constants of a biaxial sample 
and of an uniaxial sample under suitablly chosen geome- 
try and boundary conditions, in which the biaxial elastic 
constants should be comparable in magnitude to their un- 
axial counterparts. Thus, we estimate that U ~ O.leV. 
Consequently, we see that a stable polar biaxial phase 
with local C^-u-symmetry is, in principle, possible with 
appropriate surface stabilization. 



IV. HYDRODYNAMIC EQUATIONS IN THE 
POLARIZED PHASE 



The hydrodynamic equations for uniaxial liquid crys- 
tals have been rigorously derived, on the basis of the 
Poisson-bracket formalism in both the nematic phase 
(Dooh) and the polarized (Coot,) phase [3(||43j]. It has 
been shown that the dynamical equation for the director 
field retains the same form in both phases, and the effect 
of polarity manifests itself only in a different microscopic 
expression of the flow- alignment parameter, A [36|. For 
systems exhibiting macroscopic phases of biaxial symme- 
try, the identification of all three spatial directions (ver- 
sus just one in the unaxial case) leads to three distinct 
flow-alignment parameters. Any system that breaks the 
same continuous symmetry should obey hydrodynamic 
equations with identical structures [56| . Thus, contin- 
uous symmetry of a system dictates the nature of its 
hydrodynamic variables as well as the characterization 
of propagating and diffusive modes, in contrast to dis- 
crete symmetries which determine issues such as the cou- 
pling of modes. In the isotropic phase with full rotational 
and translational invariancc, if we restrict our attention 
to isothermal processes, the hydrodynamic fields consist 
only of the four conserved variables: number density and 
the three components the momentum density (we do not 
consider energy). In the biaxial phase, we have three 
additional hydrodynamic fields due to broken rotational 
symmetry. Therefore, the total set of hydrodynamic field 
variables, {$ a (r, i)}, in the biaxial case, consists of the 
number density p(r), momentum density g(r), nematic 
director field n(r), and polar director field p(r). We note 
that the director fields n(r) and p(r) are not independent 
but satisfy the orthogonality constraint n ■ p = 0. The 
conserved and broken-symemtry variables evolve in time 
according to [H, lilj 



a$„(r,t) 
dt 



v a *• at 



SF 
M> h (r) 



(26) 



where T a b is the symmetric dissipative tensor, F is the 
total free energy shown in Eqs (f23|) and ([24]) . and V a is 
the reactive term or nondissipative velocity given by 



V a (r) = - J d 3 r' {$ a (r),$ b (r')} 



dU 



(27) 



where the Poisson bracket between the field variables in 
Eq. ([27]) is defined as 

5$ a (r) 0$ 6 (r') 



d 



an 



Op 



OifJ, 



(28) 



and can be computed as coarse-grained expressions of 
suitably defined microscopic variables. In the following, 
we will present a summary of the necessary ingredients in 
the derivation of hydrodynamics via the Poisson-bracket 
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method. The detailed evaluation of the various Poisson 
brackets can be found in the Appendix. 



Poisson-bracket Formalism 



Using Eq. (|28[). it is straightforward to evaluate the 
Poisson-bracket relations between the various hydrody- 
namic fields. Our goal is to obtain equations describing 
the dynamics at long wavelengths. We will therefore only 
keep terms of lower order in the gradients of the hydro- 
dynamic fields and expand the 5-function as 



5(r 



6(r 



- Ar^V fc 5(r 



-0(V 2 ) 



(29) 



where Ar QM = r QM — r". In addition to the microscopic 
definitions of polarization and alignment tensor in Eqs. 
([7]) and (TTTj) , we have the microscopic mass density, p(r) 
and momentum density g(r), defined in the usual way 



p(r) = ^m^r-r^) , 
g(r) = ^p a ^(r-r Q ^) 



(30) 
(31) 



where p QM is the momentum of the /i-th atom on the 
a-th molecule. The macroscopic mean variables describ- 
ing the dynamics of equilibrium fluctuations in the sys- 
tem are obtained from the microscopic ones after coarse- 
graining, p(r) = [p(r)] c , g(r) = [g(r)] c , as described, 
for instance, in Ref. [431 ]. The required non- vanishing 
Poisson-brackets are: 



{p(r),9i(r')} 
{Qij(r),g k (r')} 



Vi6(r - r')p(r') , 
5i (r')V^(r-r') 
-V< [8(r-v') 9j (v')] 
tV*Q«(r)]*(r-r / ) 
-Ay H (r)V^(r-r') 



(32) 
(33) 
(34) 



where 



SikQji + SjkQu 
I 



-z^ijQki 



T^SklQi 



6AI 

+0{Q%) , 



SikSji + SuSjk — -SijSki 



(35) 



in which we average both the trace and anisotropy over 
the system and thus sum over the index a in Eqs. (|19p 
and (|2"0|) . The Poisson-bracket relation shown in Eqs. 
([34]) and ([35]) has been previously obtained by Stark and 
Lubensky in [43| for uniaxial nematic liquid crystals. In 
their expression of Xijki, they also obtained a second- 
order term in Qij by the inclusion of the dynamics of the 



amplitude variable of Rij , which ultimately did not con- 
tribute to the hydrodynamical equation for the nematic 
director. For biaxial liquid crystals, both Eqs. ([3~4"|) and 
(|35[) retain the same form when we use our molecular 
moment-of-inertia expression [Eq. (|15[) ] in the Poisson- 
bracket derivation. 

The director fields themselves have no microscopic def- 
inition in terms of the canonical coordinates and mo- 
menta of the individual atoms. Their Poisson-brackets 
must therefore be obtained indirectly from those of the 
alignment tensor Qij. Using the definition of the polar 
and nonpolar directors given in Eq. ([1]) and the chain rule 
of derivatives we obtain, respectively, 



1 

2T 



S + T 



S-T 
x {Qki(r),gj(r')}pi(r) 



ni(r)n k (r) 

(36) 



{m(r),gj(r')} = - 



*£"(«•) + 



T 



S-T 



Pi(r)pk(r) 



where 



x {Qkl(r),9j(r')}ni(r) 



$ij( r ) = Sij(r) - pi(r)pj(r) , 
Sff(r) = 6ij(r) - ni{r)rij{r) 



(37) 

(38) 
(39) 



are the projection operators onto components perpendic- 
ular to p and n, respectively. Upon insertion of Eq. (|34[) . 
we obtain 



{Pi(r),9j(r')} 



[V jPi(r)} 5 (r-r' 
xVJ(r-r') , 



2^ ( S fk( r ) - ^niWn*Wj ^kijmPi(r) 

(40) 



{ ni (r), 9 j(r')} = [V jPi (r)} 5 (r-r') 



(r) 



Tpi(r)pk{r) Tmi(r)m k {r) 



S-T S+T 
x\ kljm ni(r)\7 m 5 {r - r') . 

B. Nondissipative Velocities 

For biaxial liquid crystals, the nondissipative velocities 
associated with the polar and non-polar director fields 
follow directly from Eqs. |36|) and (|37l) and takes the 
following respective generalized forms [23, [2!^ [3l[ : 

K P = -VjVjPi + \ p ijk V k Vj , (42) 
C = -VjWjUi + \? jk V k Vj , (43) 



(41) 
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where 



A ijk 



- [(Ai + l)p 3 m k 4 (Ai - l)p k mj] rrii 
-\\{^ + l)Pjn k + (A 3 - l)p k nj]ni , (44) 
i [(A 2 + 1) rrijUk 4 (A 2 - 1) m k nj} nii 

+ l)Pjn fc + (A 3 - l)pknj]pi • (45) 



In Eq. (J44J) , only one independent reactive coefficient 
that is perpendicular to both n and p enters into the 
nondissipative velocity of the polar director field p due to 
the orthogonality constraint p • n = 0. As shown in the 
Appendix, these flow-alignment parameters Ai, A2 and 
A3 can be related to the order parameters and molecular 
parameters I and AI as follows: 



Ai 
A 2 
A3 



J_J_ 1 S_ 

&TAJ ~ 2 ~ 3T ' 
I I 



3(5 + T) AI 



1 



S 



I 



3(5 - T) AI 



- 1 



3(5 + T) ' 
5 

3(5 — T) ' 



(46) 
(47) 
(48) 



In the limit of T = where the m and p directions are no 
longer distinct and where the symmetry of our mesogenic 
molecules becomes uniaxial, the flow- alignment parame- 
ter Ai drops out of our model, and both A 2 and A3 reduce 
to the expected uniaxial value °fA = -|(l4^y^) [43l |. 

For the conserved hydrodynamic variables, the mass- 
conservation law implies for the nondissipative velocity 
of mass density 



V" 



(49) 



Incorporating the distinction of the three spatial direc- 
tions, the non-dissipative velocity for the momentum 
density takes the following form: 



9i9j 

i 

P 



+V r , 

+Vr, 



„ SF rrT . SF rrT . SF 

Sp Spj Srij 

S + T \ SF~ 

TpjPk _ Tmjm k \ 5F_ 
S-T S + T J^^Snj 



(50) 



C. Dissipative Terms 



Now we turn our attention to the dissipative terms. 
The tensor of viscosities will be more complicated due 
to biaxial symmetry leading to nine independent viscosi- 
ties (33j . In principle, the lack of p — » — p symme- 
try would permit new terms to appear in the viscos- 
ity tensor. However, these new terms appear only in 



terms higher-order in gradients. Thus, just as in the 
uniaxial case, the viscosity tensor for the polar biaxial 
phase would be the same, to the lowest order, as that 
for the non-polar counterpart with local D 2 ^-symmetry. 
Since the form of the viscosity tensor rjijki , for which 
a ij = VijkiU k i, where a\j is the dissipative stress tensor 
and U k i = h (V k vi + V/Ufc) is the strain-rate tensor, for 
the biaxial nematics is well known, we shall simply quote 
the result from literature [3ll (331 : 



Vijkl 



aiPiPjPkPi 4 a 2 n i n J n k n l 4 a 3 <5™(^ 
+a 4 {PiP n k ni +p k pin l rij) 
+a 5 (piPjSft + PkPiSfi) 
+a 6 (ninjSu + n k rnS%) 

+a 7 (pjPiriin k +p ] p k n i n l + PiPkUjUi 4 PiPin 2 n k ) 
+a 8 (pjpi6% + p 3 p k 5% + PiPkSJi 4 PiPid$) 
+a 9 (njniS" k + njn k 5™ + riin k S"l 4 runiSjl) . 

(51) 



where 



$ij = - PiPj ~ n * n 3 ( 52 ) 

is the operator projecting onto components along m 
p x n. Therefore, 1 
the following form: 



p x n. Therefore, the dissipative stress tensor takes 



aipiPjp k piU k i 
-azS™p k piUu - 
+a i p k pin i n j U k i 
+a 6 n i n j S k n l U k i 4 
+a 7 p J p k n l n l U k i 
+a 8 pjPiS^U k i 4 
+a$p i piSJ l k U k i + 
+a 9 rnn k SJ}U k i 4 



a2niTijn k niU k i 
asSTjn k niU k i 4 
4- asPiPjSuUki 



h a 3 5%U kk 
aiPiPjn k niU k i 



abPkPiS^Uu 
a 6 n k rnS™U k i 4 a 7 p 3 pinin k U k i 
+- a 7 p l p k n j niU k i 4 a 7 piPin 3 n k U k i 
asPjPkSaUki 4 a 8 p t p k SJ}U k i 
agrijUiS^Uki + agrijUkS^Uki 
agmriiSJlUki . 

(53) 



For the polar and nematic director fields p(r) and n(r), 
we must include all dissipative terms allowed by the sym- 
metry. The odd symmetry of dissipation under time- 
reversal implies that only dissipative couplings to SF/Spj 
and SF/Srij are allowed. In view of Eq. (|2"6"|) . the dissi- 
pative part of the polar- and nematic-director equations 
can be schematically written as follows: 



(54) 
(55) 



dpi _ 


SF 


SF 


dt 






dnP _ 


r „n SF. - 


SF 


dt 




ij s Pj 


coefficients 







\ L ij J ' 

are the rotational viscosities. The Onsager prin- 



(iT) 

ciple requires another constraint that 



ij ' 



Im- 



posing the unity and orthonormality conditions on the 
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director fields, n 2 = p 2 = 1 and n ■ p = 0, we obtain a 
total of twelve independent rotational viscosities for the 
dissipative tensor: 



dt 



-rrii 



HP] 



SF 
8n« 



SF 

, SF 
Snj 
SF 



5F 

72^7 1- 

'^Sn~- 

SF 
■Isrij-z— 4 



SF 
, SF 



Is 7 



SF' 
lerrij — 
SPj. 



(56) 



~J ( n kPj + n jPk) U kj ni 



IP] 



SF 

Srij 



SF 



SF 



SF 



SF 
SF 



'IP] j Ifij- \-jrrij — 

5 Pj Sp-j S Pj 



SF 



SF 



SF 



Sn^ 
SF 



SF 



S Pj 5 Pj 



5n 
5F_ 
j Sp~ 



(60) 



drf 
dt 



-HP] 



SF 
Srii 



SF 

8ri4 



73 mj 



SF 
Srii 



liP] 



SF_ 

Sp 3 
SF 



5F_ _ 

' Sp-j 
SF 



73"lj 



SF 
Spj 
SF 



-rrii 



liPj- 



-I2P] 



5F_ 

Spj 



7 x n 



Snj 
SF 



Hp, 



'73 r 



5n 3 

SF 
>Sp~ 



(57) 



We can easily see that Eqs. (j56| and (|57|) readily satisfy 
the constraints that (drii /dt) — pi (dpi /dt) = and 
Tit (dpi/dt) + pi (drii/dt) = 0. For simplicity, we will set 
all twelve rotational viscosities to be 7 in our subsequent 
calculations. As in the case of polar unaxial liquid crys- 
tals, the only difference between the two sets of equations 
come from differences in the free energy. 

We now collect the reactive and dissipative terms and 
write down the hydrodynamic equations for the field 
variables. We will then consider linearized fluctuations 
about the ground state of uniform director fields along 
the n and p-directions as well as of a uniform density 
Po. Inertial terms connected with rotations of align- 
ment axes are generally negligible for unaxial systems, 
and the same will be true for biaxial nematics. The in- 
ertial terms connected with flow should be the same as 
those for normal liquids, and they can, of course, not 
be generally neglected even though in many cases they 
will be small against friction [13]. Defining the vorticity 



Ujj 



h (VjUj — Vji>i) in the director equations, the full 



hydrodynamic equations are as follows: 



dp 
at 



-V-g(r) , 



(58) 



dt 



„ 9i9i „ SF ! SF r _ . SF 

■Vj— - pVi — + [ViPj] j— + [ViTij] j— 



Sp 



S Pj 



, Tp S + T \ SF 

OA - — —TljUk AklimPlT 

J $P] 



1 



2T\i k S-T 



TpjPk Tmjm k \ SF 



1 



S \ jk S — T S + T 



]"lj ) 



' Srii 



(59) 



-—- = -VjVjUi - uoijUj + (m k n 3 + m,jn k )U k jm. 



A, 



dt 



+Y (P* n i +P] n k) U k3 Pi 



-Pi 
-Pi 



-IPi 



IP] 



SF 
Snj 
SF 
5 Pj 
SF 
Snj 
SF 



SF 



JTlj- 



SF 
Srij 



SF_ 

SF 



^sfr^wr^spi 



5F_ 

SF 



Srii 



■1 m j 



Srij 



SF 



Sp, Spj 



where 



SF 
5p 3 



•frrij 

SF ... (Sp\ B 



SF 
Spj 



Sp \PoJ Po 



D 

Po 



(61) 



(62) 



Vjp - K pl Wj(V ■ p) + K p2 {V x (V x p)). 



SF 

Snj 



-(K p2 - K p3 )p s (V s pk) (Vj-pfc) 
+ (K p2 -K p3 )V b (p bPs (V s p 3 )) + 0(V 3 ) , (63) 

= -K nl Vj{V • n) + K n2 {V x (V x n))j 

-(K n2 - K n3 )n s (V s n k ) [V n k ) 

+(K n2 - K n3 )V b (n b n s (V s ^)) + O (V 3 ) . 

(64) 

The omitted terms in the expressions for SF/Sp, 5F/5pj, 
and SF/Snj do not contribute to the linear regime of in- 
terest. In spite of the apparent complexity of the biaxial 
elastic free energy in Eqs. (|23f and I] 24 p. most of the 
terms are highly nonlinear and, fortunately, do not come 
into consideration in the linear hydrodynamics regime. 



V. HYDRODYNAMIC MODES IN THE 
POLARIZED PHASE 

In the isotropic phase with full rotational and transla- 
tional invariance, if we restrict our attention to isother- 
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mal processes, the hydrodynamic fields consist only of the 
four conserved variables: number density and the three 
components the momentum density. The corresponding 
four hydrodynamic modes are the two propagating sound 
waves associated with the decay of density and longitudi- 
nal momentum fluctuations, as well as two shear modes 
describing the decay of the two transverse components 
of the momentum density. In the less-symmetric unaxial 
nematic phase, be it macroscopically polar or otherwise, 
there exists two additional modes associated with the 
transverse fluctuation of the director field, which arises 
from the residual unbroken S f O(2)-symmetry. The sys- 
tem now possesses six hydrodynamic fields in total, and 
the modes now separate into a set of two diffusive modes, 
associated with the velocity and director fluctuations 
transverse to both the wavevector and equilibrium di- 
rector orientation, that characterize the coupled decay of 
vorticity and twist fluctuations [36|, [57j , and a set of four 
remaining longitudinal modes that consist of two prop- 
agating sound waves and two diffusive modes describing 
the coupled decay of splay and velocity, longitudinal to 
the wavevector, fluctuations [57| . In a nematic system 



that exhibits macroscopic polarity, the sound waves re- 
ceive corrections from the parity-odd terms in the free en- 
ergy, but the overall structure of the longitudinal modes 
remains qualitatively the same as the non-polar phase. 

We now consider the long wavelength, low frequency 
dynamics in a uniformly polarized biaxial phase, char- 
acterized by uniform equilibrium values po of the num- 
ber density, vo = of the flow velocity, po = y of 
the polar director field, and no = z of the nematic 
director field. In order to evaluate the hydrodynamic 
modes of the polarized biaxial state, we expand the hy- 
drodynamic equations to linear order in the fluctuations 
of density, Sp = p — po, momentum Sg = pov, polar 
director field 5p = p — y, and nematic director field 
Sn = n — z about their equilibrium values. Since p and 
n are orthogonal, there are only three independent vari- 
ables associated with the broken symmetries, chosen as 
Sn x , Sn y , and Sp x (Sp z = —Sn y ). Again, we examine 
only isothermal processes. Considering wave-like fluctu- 
ations whose space and time dependence is of the form 
exp [iq ■ r — iuxt] , we obtain the linearized equations for 
the polar biaxial phase: 



dtSpd = -ipo q ■ v q , 



(65) 



d t 8p x = Sp q 



* — (Qx + q y + q z 

Po 



+ 5pl 7 \-K& (q)q 2 + AtfW (q x q z + q x q y ) 



+<5< 7 [-K x n) (q)q 2 + AK™ (q x q z + q x q y ) 
+Sn«j \K^(q)q 2 - K y n \q)q 2 - (q x q z + q y q z ) + AA» (q x q y + q y q z ) 



1 + Ai 



1- Ax 



(66) 



d t 5n x = Sp q 



i — [Qx -q y - qz) 

Po 



-K<?\q)q 2 - AK^ (q x q z +q x q y ) 



+<5n q7 \-Ki n \q)q 2 + (q x q z + q x q v ) 



+5n v 7 



K^\q)q 2 - K y n \q)q 2 - AK^ (q x q z - q y q z ) + AK^ (q x q y + q y q z ) 

1-A 2 



1 + A 2 
' i — - — i q 



qx 



(67) 



d t 5nl 



$Pq 



i — {q z + q y - q x ) 
Po 



+ dpi 7 \K<f\q)q 2 + AifW (q x q z + q x q y ) 



-8f%1 [K x n \q)q 2 - AA» (q x q z + q x q y ) 

-<5n q7 \K^(q)q 2 + K y n) {q)q 2 + AK^ (q x q z - q y q z ) - AA» (q x q y + q y q z 



r./i + A 3 \ 






K 2 ) q \ 


-< 


K 2 H 



(68) 
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Pod t vl 



+Sp^ 



+c5n q 



■ Ci 2 , / 1 + Ai 

-i — q x - iC 2 p q x q + — ~ — 
Po V 2 



q x q y 



1 + A 2 



+5n q 



Ki n) (q)q 2 q z 



1 + Ai 



, ; J AK^q x q y q z 



q x q y q z 



Bq x q z + z 

["39a + "8^ + aggf] - [("5 + «s) q x q y ] - Wq [(a 6 + ag) fefe] , 



1 + A 



— ) AJfW 

2 y 



(69) 



PO^q = (5/?q 

+ ^q 

+<5< 



+<5n q 



■ Ci 2 , f l-A 3 

-t — g a - iC 2 poq y q + — - — 
Po V 2 



1-Ai 



B 



Bqxqv-il?-^)^®^-, 



Po/ V 2 / \po 
1-A 3 



2 / 



AK^q x q 2 z 



Bq y q z 



1 " M AK^qlq z + i( l -±^)K^{^q z 



2 y z ^ ™ \ 2 / m \ 2 

Wq [(a 5 + a 8 ) - w q [ai^ + a 7 g£ + a 8 <? 2 ] - « q [(«4 + & 7 ) q y q z ] , 



(70) 



Pod t Vq 



+s P x cl 



-1 — q z 

PO 
Bq x q z 



iC 2 p q z q 2 
I + A3 



1 + A, 



q y q z 



j AK^q x q y q z 



1 - A 5 



K x n \q)q 2 q x +i 



1 - A, 



q x q y 



IJ q - + , ( I±*») ^(q)^ + i (i^) ARM 



qlq y 



-v* [(a 6 + ag) - [(a 4 + a 7 ) q y q z ] - w q [a 2 q z + a 9 q 2 x + a 7 ^] 

I 



1-Af 



K v n) m 2 q y 



(71) 



where we have defined the following quantities: 

JfW(q) - [/fm^ + ^^ + ifnS^/g 2 , (72) 

Jf W(q) = [Jf„i<£ + K n2 q 2 x + X„ 3 <Z 2 ] /<? 2 , (73) 



^ P) (q) 

^)(q) 



[i^ P l9 2 + K p2 q 2 z + K p3 q 2 y ] /q 2 , (74) 
[K pl q 2 z +K p2 q 2 x + K p3q 2 y ]/q 2 , (75) 



AA"(P> = K p2 - K pl , 
ARM = if„2-^„i 



(76) 
(77) 



The system has seven hydrodynamic modes described 
by dispersion relations «j(q). Stable modes describ- 
ing decaying fluctuations correspond to 3?[iw(q)] > 0, 



while a negative sign in 3?[zw(q)] signals the instabil- 
ity of the corresponding uniform state. The seven hy- 
drodynamic modes are all coupled together. There 
is no longer a separation between sets of transverse 
and longitudinal modes. In terms of the basis of field 
variables (Spq, <5p q , <5n q , <5n q , u q , v^, v q ), probing their 
mode structures amount to diagonalizing a linear sys- 
tem of seven coupled differential equations. While the 
numerically obtained expressions of the normal modes 
and their associated eigenfrequencies, expanded to 0(q 2 ), 
are rather complicated, the salient feature of the normal 
modes in the polar biaxial phase is that the new parity- 
odd terms, proportional to the elastic constant B, give 
quantitative corrections to the sound modes as well as to 
the diffusive modes. In total, there are a pair of sound 
modes with the usual (when compared with other fluids) 
sound speed of i^/Ci/po and a complicated anisotropic 
decay rate, and five modes with dispersion of order 0(q 2 ). 
Comparing with the hydrodynamic modes found for un- 
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axial liquid crystals, there always exists a new purely 
diffusive mode which may be identified as the Goldstone 
mode of the uniaxial-biaxial transition [33|. Overall, the 
diffusion peaks contain angular dependence as all three 
spatial axes are now identified in the polarized biaxial 
phase. In particular, we will express our results in terms 
of spherical coordinates: q x = sin 8 cos 0, q y = sin 8 sin 0, 
and q z — cos 8. Expanding the hydrodynamic frequencies 
us to quadratic order (icj = rjiq + i]2q 2 ), and solving the where 



eigenvalue problem perturbatively, we obtain the follow- 
ing explicit expressions for the propagating sound modes: 



±iq\ —-q 2 -D(<l>,9) 



(78) 



D(0,0) 



1 B 2 
8Ci 



[-1 - 3 008 26* + 2 sin 2 8 (cos 20 + sin 20) + 2 811126* (2 cos - sine 



+i>2 cos 4 8 + sin 4 8 [z^ cos 4 + v\ sin 4 + 2{y§ + 2i/g) cos 2 sin 2 0] 

+ cos 2 8 sin 2 8 [2(v e + 2v 9 ) cos 2 + (z/ 4 + 3vj) sin 2 + (z/ 4 + v 7 ) cos sin < 



(79) 



where we have defined the kinematic viscosities as V\ — 
otijlpQ as they are the natural quantities that deter- 
mine momentum relaxation. In particular, along the 
three specified directions along q x where (0, 8) — (0, 7r/2), 
along q y where (0, 9) — (7r/2, 7r/2), and along q z where 
8 = 0, the diffusive components of the first sound mode 
reduce to the following limits: 

p(0 = M =|) = (g! + , 3 )^ (80) 
v(j>=^G=l) =v x ql (81) 
I? (0,0 = 0)= Lg-+^g 2 (82) 

The sound modes are anisotropic, as expected from the 
calculation of the structure factor and the free-energy 
analysis in section II. In addition, the coupling of density 
to splay fluctuation can render the sound modes unsta- 
ble for certain direction of propagation, although viscous 
dissipation tends to stabilize them. Assuming, for sim- 
plicity, that Vi — v for all values of i, we find that the 
fastest growth of fluctuations occurs for q in the zz-plane, 
corresponding to q y = or = 0. In this case we obtain 

D (0 = 0, 9) /v = 1 + sin 2 29 + a( sin 29 - cos 29) , (83) 

where a = ^B j (2vC\). If T*(<f>,9) < the modes are 
unstable. This occurs provided a > 1. For small values 
of a — 1 the sound waves are unstable for all angles 9 < 
9 c (a), with 9 C w (l + a)/2. 

The remaining five modes all have dispersion relations 
that are at least of order 0(q 2 ). The corresponding nor- 
mal frequencies can be obtained from the solutions of 
a fifth-order equation. Their exact expressions are quite 
complicated, and we will specialize to a few limiting cases 



below whose physical properties can be readily illumi- 
nated. 

In the incompressible limit, the sounds modes will have 
leading contribution from 0(q 2 ) just as what we found 
for the polar unaxial case, as terms proportional to B 
(and D to higher orders) vanish. The vanishing of these 
partiy-odd terms can be readily seen in the linearized 
regime in Eq. (|79[) of which the q 2 terms contain contri- 
butions quadratic in B. The physical nature of this lead- 
ing O ((^-contribution will depend on the relative mag- 
nitude of the various viscous coefficients. To illustrate 
this point concretely, we now consider the one-constant 
approximation, often employed in numerical simulations 
of liquid-crystal systems, in which the corresponding 
mode structure simplifies considerably. In this approxi- 
mation, we neglect the anisotropy in the elastic constants 
along both polar (p) and non-polar (n) directions, setting 
K p i = K p2 = K p3 = K p and K nl = K n2 = K n3 = K n . 
We will consider, in the subsequent subsections, firstly 
the three limiting cases in which we set K p = while 
keeping a nonvanishing K n with linearized fluctuations 
along only q = q x , q = q y or q = q z , separately. We 
will then consider three additional limiting cases, where 
fluctuations are again only found along each of the three 
distinct directions but with K n = and a nonvanishing 



A. K p = 

When we set K p = 0, the mode structure becomes es- 
pecially simple when we consider only fluctuations along 
q y . In this limit, we obtain the following expression for 
the quintic algebraic equation that governs the 0(q 2 )- 
contributions to the eigenmodes (iu> = rjiq + r]2q 2 ): 
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m (m - 2^ 8 ) 



dril + 2v 7 C x ril + 



CiK n 
4p 



(i-A 3 ) \ 2 m + 



r 



1C1KI 

4p 



(1 - V 



(84) 



As expressed in Eq. (|84|) , we readily see that there are 
two modes that decouple from the rest. The first is 
the diffusive mode described by the dispersion relation 
iui = 2q 2 i/g that controls the decay of v x , which decou- 
ples in this limit. The second mode has dispersion rela- 
tion iu = and describes the relaxation of Sp Xl which 
does not decay since we set K p = 0. The remaining 
three modes are governed by the solution of the cubic 
equation in square brackets in Eq. (|54"|) . Using Math- 
ematica, we see that they consist of a pair of complex 
conjugate solutions corresponding to propagating modes 
with speed of propagation of 0(q 2 ). The third solution 
is real, corresponding to a diffusive mode. The nature of 
these solutions, and hence the stability of the three cor- 
responding modes can be examined using the Descartes' 
rule of signs from algebraic theory: since all four coeffi- 
cients in the cubic equation have the same positive sign 
(1/7 > 0), there are no positive real zeroes to the cu- 
bic equation. Likewise, when we substitute 772 — > — 772 
in the cubic equation, we observe three sign-changes be- 
tween the four terms, when written in descending alge- 
braic power, and the Descartes' rule of signs indicates 
that there would consequently be either three or one real 
negative roots to the cubic equation. Since we already 
know that two of the three roots are in general complex 
conjugates, the remaining one real root must hence in 
general be negative and, thus, corresponds to a physi- 



cally stable eigenmode. The stability of the propagating 
modes, on the other hand, is determined by the proper 
sign of their real part, and an analysis of its precise na- 
ture can only be preformed with the knowledge of its full 
expression. 

The mode structure for uni-directional fluctuations 
along the other two directions, along q = q x and q = q z , 
are both considerably more complicated, as the equa- 
tions governing the 0(g 2 )-contributions to their respec- 
tive eigenmodes are full quintic equations. Moreover, the 
nature of the algebraic roots and the corresponding sta- 
bility of the eigenmodes would, in general, depend upon 
the relative magnitudes of the various elastic and vis- 
cosity coefficients. The detailed elucidation of the gen- 
eral hydrodynamic mode structure of the biaxial phases, 
formed by bent-core nematogens, will be presented else- 
where [58j |. 



B. 



K„ = 



Analogously, we can consider the simplified limit of 
nonvanishing K p and set K n = 0. In the case of fluctua- 
tions occuring only along q = q Xl we also find a compli- 
cated mode structure governed by the following factor- 
izcd quintic equation: 



J 



im - 2zaj) { - C^l + (B 2 1 - 2i*Ci) 4 




4/90 

(1-Ai) 2 

4/90 



{B 2 -C l K p )-2 1 B 2 u s m 



(B 2 - C\K p ) K pl \=Q 



(85) 



In this limit, v z is decoupled, and its decay is described 
by a diffusive mode iu> — 2q 2 vg. There is a mode iu> = 
corresponding to the decay of Sn x , which does not de- 
cay since we set K n = 0. Of the remaining three modes 
given by the solutions of the cubic equation contained in 
the curly brackets of Eq. |85|) . two are complex conju- 
gates and one is real. The stability of the real mode is 
controlled by the relative magnitude of the various vis- 
cosity and elastic coefficients, one of which includes the 
factor of {B 2 — K p C\). Comparing our free energy anal- 
ysis in Section II and the examination of the eigenmode 
equation in Eq. (|85|) , we find that the condition signal- 



ing the onset of splay instability as found in Section II, 
namely, \B\ > w 1 K p C\ may actually be a necessary but 
not sufficient requirement for mode instability. Its rig- 
orous confirmation would require the full analysis of the 
eigenvalue problem [58|. 

The mode structures for the other two limits, in which 
fluctuations are found along only the directions of q = q y 
or q = q Zl are again governed by the full quintic algebraic 
equation. Needless to say, the mode expressions become 
even more cumbersome when we re-introduce anisotropy 
in the elastic constants, though they all share several 
overall features. In general, the precise physical nature 
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of all five modes, whether diffusive or propagating, de- 
pends on the relative magnitude of the various viscosity 
coefficients and elastic constants. When the solutions to 
the full quintic equation, and thus the corresponding hy- 
drodynamic modes, are complex conjugates, they signify 
that the leading contribution to their propagating com- 
ponents are of order 0(q 2 ). 
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APPENDIX A: DERIVATION OF 
FLOW-ALIGNMENT PARAMETERS 



VI. CONCLUSION 

In conclusion, we have examined in this work the hy- 
drodynamics of a biaxial nematic system which is macro- 
scopically polar, with local C2„-symmetry. Our work is 
motivated by the recent experimental discovery of the 
long-sought biaxial nematic phase composed of micro- 
scopically polar bent-core mesogens. Our present results 
are extension to those previously obtained for the hydro- 
dynamics of polar uniaxial liquid crystals. In spite of the 
much more complicated structure of the hydrodynamic 
equations due to the complete breaking of rotational in- 
variance and to the discrete orthorhombic symmetry, we 
have found that polarity and nematic biaxiality can coex- 
ist, although the system exhibits an intrinsic bulk splay 
instability for a range of parameters and along certain di- 
rections of the wavevector. Moreover, we have found that 
the newly allowed parity-odd terms in the free energy as 
well as the orthonormality constraints between the direc- 
tors in the biaxial system provide additional corrections 
to the hydrodynamic mode expressions and render them 
highly anisotropic. Polarity does not, however, alter the 
linear hydrodynamics of the biaxial phase in the incom- 
pressible limit, just as in the more symmetrical uniax- 
ial state. In deriving the hydrodynamic equations using 
the Poisson-bracket formalism, we have also computed 
the flow-alignment parameters along the three symme- 
try axes of the polar phase in terms of microscopic pa- 
rameters associated with the molecular geometry of the 
constituent biaxial mesogens. 

It is our hope that the results presented in this paper 
will contribute to the theoretical understanding of biaxial 
nematic phases, which have continually generated great 
interest from the scientific community due to their many 
fascinating properties as well as their immense poten- 
tial for technological applications. In particular, there 
has been a recent focus on the phenomenon of flcxoclcc- 
tricity (59|-an effect uniquely enhanced by the particu- 
lar boomerang shape of these liquid crystals-that car- 
ries the potential of serving as a basis for environmen- 
tally friendly micro-power generators. Furthermore, our 
study of the hydrodynamic modes would complement 
the experimental study of the dynamics of the phases 
of these mesogenic systems using techniques such as dy- 
namic light scattering [6(j. It is our further hope that 
our results, presented here in the context of equilibrium 
physics, will also serve as a starting point in understand- 
ing the behavior of nonequilibrium systems where polar- 
ity arises from actively driven mechanisms. 



In this Appendix, we present details on the derivation 
of our flow- alignment parameter expressions found in Eqs 
flU, (g7]), and flU]), defined via the following: 



1 



[(Ai + l)pjm k + (Ai - l)p k mj] m l 



-0KA3 + l)Pjn k + (A 3 - l)p k nj]m , (AI) 



A«fc = 7: [(A 2 + 1) rrijUk + (A 2 - 1) m k nj] rrii 



+ 2 [(A3 + + (A3 - l)Pkn 3 ] Pl , (A2) 



and 



------ 

1 ~~ 6T A? ~ 2 ~ 3T ' 



A 2 = 
A3 = 



1 



3(5 + T) AI 
1 I 



- 1 



- 1 



S 



3(5 + T) 
S 



3(5 — T) AI 3(5 — T) 

Starting from Eqs. (|35J) and f3"7|). 



(A3) 
(A4) 
(A5) 



1 



{Pi(r), 9j (r')} = — 



S + T 



Wr),Sj(r')} = g 



S-T 
x{Q kl (r), gj (r')} Pl (r) , 
1 



+ —PiPk 



T 



' ' S-T rit "- S + T 
x {Qki(r),9j(r')}ni(r) , 



dU 



and applying Eq. (|27j) . 

V a (r) = - J d 3 r' {$ a (r),$ h (r')}- 



9$ fc (r') 

to the director fields p(r) and n(r), we obtain 



(A6) 

TOjTOfe 

(A7) 
(A8) 



V p 



V n 



2T S-T 



rriimk , PiPk 



S+T S-T 



(A9) 

(V m 1>;) . 
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Using Eq. (|A11[) and the biaxial form of the alignment 
tensor in Eq. (jA"l~2)) 

2 2 
Ayfci = SikQji + SjkQu — -SijQki — -^SkiQij 



+0{Q%) , 



Q«W = S(r) 



Ui(r)nj(r) - -<5y 



+T(r) [p,(r)pj(r) - m^rjm^r)] 



(A12) 



(All) we obtain the following, upon further simplifying: 
I 



n-inkXkjimPj 
mim k XkjimP] 

TnimkXkjlmTlj 

PiPkXkji m rij 



( I 



\6AI 
\6AI 
V6A7 3 



- g- I n lPn 



\6AI 3 



7 ^ 2 Q) 

+ 77^ « 



6A7 3 



T -f)^+(^j- T -|]'"^' 



-5 mim r 



— : T m m ni 

6A7 3 ' 



6AI +T I ' 



(A13) 
(A14) 
(A15) 
(A16) 



r 



Further decomposing Eqs. (|A13j) - (|A16j) into symmetric and (|A10|) . we obtain the desired results in Eqs. (|A3 
and antisymmetric parts and comparing with Eqs. (|A9[) (|A5|) . 
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